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This paper describes the measurement of collective flow for 
charged particles in Au+Au collisions at ^/s^Tiv = 130 GeV 
using the PHOBOS detector at the Relativistic Heavy Ion 
Collider (RHIC). An azimuthal anisotropy is observed in the 
charged particle hit distribution in the PHOBOS multiplic- 
ity detector. This anisotropy is presented over a wide range 
of pseudorapidity (77) for the first time at this energy. The 
size of the anisotropy (va) is thought to probe the degree of 
equilibration achieved in these collisions. The result here, 
averaged over momenta and particle species, is observed to 
reach 7% for peripheral collisions at mid-rapidity, falling off 
with centrality and increasing \r)\. Data are presented as a 
function of centrality for \r]\ < 1.0 and as a function of 77, 
averaged over centrality, in the angular region -5.0< rj <5.3. 
These results call into question the common assumption of 
longitudinal boost invariance over a large region of rapidity 
in RHIC collisions. 

PACS numbers: 25.75.-q 

The study of collective flow in non-central ultra- 
relativistic heavy ion collisions is important because it 
can provide information on the initial spatial anisotropy 
of the reaction zone, conditions present in early stages 
of the collision |l]-^| and the degree of thermalization 
attained during the evolution of the collision |||| . In ad- 
dition, collective flow can affect other measurables of in- 
terest, such as two particle correlation functions and the 
slopes of transverse momentum distributions (l],Q|. Ef- 
fects of collective flow have been observed in nuclear col- 
lisions over a wide range of collision energies and species 

Hydrodynamic models, which assume local thermal 
equilibrium at all points, are generally thought to pre- 
dict maximal flow. Such models are fairly successful at 
reproducing the mid-rapidity flow results at RHIC for 



the more central events and lower transverse momenta 
[jl7],[l8]]. This implies substantial early equilibration in 
these collisions. Hydrodynamic models predict roughly 
constant flow over a broad pseudorapidity region either 
through a full three-dimensional calculation [ fl8[ or by 
assuming longitudinal boost invariance. This paper pro- 
vides an additional constraint on such models. 

The analysis presented here is based on data taken be- 
tween June and September, 2000, during the first RHIC 
physics run. Results are presented for Au+Au collisions 
at \/ s nn = 130 GeV as a function of centrality over a 
narrow range of pseudorapidity (77) , while minimum bias 
results are shown over a large 77 range (-5.0< 77 <5.3). 

The PHOBOS detector employs silicon pad detectors 
to perform tracking, vertex detection and multiplicity 
measurements. Details of the setup and the layout of 
the silicon sensors can be found elsewhere |h|-|21]. For 
this running period, detector components relevant for 
this analysis included the first six layers of both sili- 
con spectrometer arms (SPECN and SPECP), the silicon 
vertex detector (VTX) , the silicon octagonal multiplicity 
detector (OCT), three annular silicon ring multiplicity 
detectors on each side of the collision point (RN, RP), 
and two sets of scintillating paddle counters (PN, PP). 

Collisions with a coincidence of two or more signals 
in each of the PN and PP counters were selected by 
the trigger. This sample corresponded to 86 ± 3% of 
the total inelastic Au+Au cross section. The centrality 
determination for the triggered events was based on a 
truncated mean of the deposited energy in the paddle 
counters. This variable is proportional to the number of 
particles hitting these counters and is monotonically re- 
lated to the number of participants, N port . More details 
on the event triggering and centrality determination can 
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be found elsewhere |22|]. Adjustments in the procedures 
used previously were made to take into account the fact 
that data used in this analysis came from an offset fidu- 
cial volume and included periods with different magnetic 
field settings. 

Monte Carlo (MC) simulations of the detector perfor- 
mance were based on the Hijing event generator |23|] and 
the GEANT 3.21 (24) simulation package, folding in the 
signal response for scintillator counters and silicon sen- 
sors. 

The anisotropy of the azimuthal distribution of 
charged particles traversing the detector formed the basis 
for this flow analysis. Uniform and symmetric acceptance 
was beneficial in terms of sensitivity to the flow signal. 
This led to the requirement that the primary collision 
vertex fall within an 8 cm fiducial region centered at -34 
cm from the nominal interaction point, in a uniform and 
symmetric part of the OCT subdetector. 

The position of the primary collision vertex was de- 
termined on an event-by-event basis by fitting for the 
optimal intersection point of the straight tracks recon- 
structed in the first six planes of each of the two spec- 
trometer arms. In addition, the vertex position along the 
beam, z v t x , was required to be consistent with the vertex 
position as determined by an independent algorithm to 
avoid pathological vertex reconstructions and reduce the 
potential for systematic effects. The second algorithm 
determined z vtx as the z position of the maximum of the 
azimuthally averaged hit density in the OCT subdetec- 
tor. Finally, the reconstructed transverse vertex posi- 
tion was required to be within 2er of the average position 
(beam orbit). 
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FIG. 1. The event distribution as a function of N 



part 



for 



triggered events (upper curve) and for data accepted for use 
in the final analysis (lower curve). 

Trigger selection yielded 1.37 x 10 6 events. A total of 
13,644 events survived the vertex cuts described above, 
yielding a global acceptance for events used in the analy- 
sis of The number of triggered and accepted events 
are shown as a function of N part in Fig. |l|. The vertex po- 
sition resolution for accepted events was estimated from 



simulated data to be ~1 mm in x, ~1 mm in y and 
~3 mm in z. 

The raw data for this analysis came in the form of 
energy depositions in individual detector pads, known 
as hits. The hit energies were adjusted for variations 
in silicon thickness and converted to dE/dx using the 
expected path length, assuming each hit to come from 
a particle emanating from the reconstructed event ver- 
tex. Pads with energy depositions greater than 0.625 of 
the peak of the minimum ionizing particle distribution 
(i.e., >50 keV) were taken to represent points of charged 
particle transit and were used in the analysis. At the 
ends of the octagon, where a single particle track often 
passes through more than one detector pad, the energy 
signals in adjoining pads consistent with being from a 
single track were added together. To avoid biases in- 
troduced by malfunctioning pads, the signals from such 
pads were replaced with signals from corresponding mir- 
ror image pads in rj, making use of the fact that the 
data were taken during symmetric collisions. This was 
done for less than 3% of pads in the OCT and less than 
1% of the pads in RN and RP. The position of each hit 
was smeared randomly with a flat distribution within its 
hardware pad boundaries and mapped into rj — <f> coordi- 
nates. The analysis was done on an event-by-event basis 
in r\ — <f> space. 

The strength of the flow is given by the n th Fourier 
coefficient of the particle azimuthal angle distribution, 



dN 



d(<f>-ih) 



1 + Y~] 2v n cos[n(cf) — ip 2 )] 



(1) 



where ip2 is the event plane angle, which is the best esti- 
mate of the azimuthal reaction plane angle, ipR, defined 
by the impact parameter and the beam axis. This anal- 
ysis was confined to n=2, the so-called elliptic flow. It 
generally followed the scheme proposed by Poskanzer and 
Voloshin ]2j| and was based on the correlation of hits in 
one part of the detector, known as a subevent, with the 
event plane angle as determined by hits in a different 
part of the detector (a different subevent). ip2 in a given 
subevent, 'a', was determined by 
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-tan 



J2i WiCOs(2<f>i) 



(2) 



where the weights, Wi, were selected to maximize reaction 
plane resolution by adjusting for acceptance and occu- 
pancy effects as described below. The sums ran over all 
hits in subevent a. 

The ip2 distribution as a function of azimuthal angle 
should be flat in the absence of detector effects. Struc- 
tures in the raw </) hit and tp2 distributions were un- 
derstood qualitatively from studies of simulated data as 
coming from gaps between sensors, channel-to-channcl 
differences in <f> phase space coverage, and the asymmetric 
production of background particles. These effects were 
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removed (and the ip 2 distribution flattened) in each 77 
annulus through the use of individual hit weights, w", 
which were proportional to the inverse of the average 
number of hits in each pad. Residual effects due to the 
variation in the vertex position were absorbed in the sys- 
tematic error. It should be noted that the final results of 
the analysis were insensitive to this weighting and that 
results with no acceptance weighting were consistent with 
the observations reported here. 
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FIG. 2. Reconstructed event plane angle correlation be- 
tween subevents OCT- and OCT+ as a function of centrality 
bin. This quantity is directly related to the event plane reso- 
lution, as shown in equation (5). 

The finite pad size in the detector led to an occupancy- 
dependent hit saturation that reduced sensitivity to flow. 
This effect was accounted for in the analysis by weighting 
the hits in a given 77 — </> bin by the average number of 
tracks per hit pad, or occupancy, calculated individually 
in different sections of the detector. The occupancy was 
determined on an event-by-event basis from the number 
of occupied (N occ ) and unoccupied (N„ nocc ) pads in each 
section. The occupancy weight in a given section was 
determined assuming a Poisson statistical distribution as 

PI 



0cc(?7, 4>) = 



li 



1 



(3) 



where /i=ln(l + N occ /N„„ occ ) is the average number of 
tracks per pad. This occupancy was used in concert with 
the acceptance weight to produce the overall weight, 



w 4 = WiOcc(r)i, (pi), 



(4) 



which was used in the determination of t/? 2 in equation 
(2). The event plane resolution, R, was calculated sep- 
arately for each centrality bin using the subevent tech- 
nique p5) as 



R = 



sj< cos[2{r 2 -^ b 2 )] >, 



(5) 



where the superscripts denote separate subevents within 
a given event and the averaging was done over all ac- 
cepted events in a given centrality bin. In this analysis, 



equal multiplicity subevents were defined by dividing the 
event into two separate rj regions, OCT-, which extended 
from -2.0 to -0.1 in 77, and OCT+, which encompassed 77 
between 0.1 and 2.0. The gap between the two angular 
ranges in the OCT subdetector was introduced to reduce 
effects due to short-range non-flow correlations between 
hits in different subevents. The final V2 determination 
was found to be robust against the choice of subevent 
used in the evaluation of the event reaction plane so long 
as the chosen subevent contained sufficient statistics that 
the reaction plane resolution could be well determined. 
In addition, the gap width between subevents was varied 
from 0.2 to 1.0 in 77 and the change in v 2 was incorpo- 
rated in the final systematic error for the analysis. The 
correlation between the event planes in the OCT- and 
OCT+ subevents is shown in Fig. ||. 




FIG. 3. Fully corrected measurement of elliptic flow, V2, as 
a function of the number of participants for 7/ <1.0. The 
black error bars represent the la statistical errors and the 
gray bands give a measure of the systematic error for each 
point at ~90% confidence level. 

The observed, resolution and occupancy corrected 
value of V2, v 2 bs , was calculated in bins of centrality and 
77 from the 77 — <p hit map using 



obs _ , < w t cos[2((j) - ip 2 )} > 
2 ~ { R 



(6) 



where the averaging in the numerator was done over the 
hits in one event and the averaging of the fraction was 
done over all events in the given centrality or 77 bin using 
the appropriate value of R for each event. 

In this analysis, vg 65 was calculated in several regions 
of the detector using the <j> position of the hits in the 
77 — (f) space. Two regions at high 77 corresponded to 
areas covered by the ring detectors, RN and RP, and 
encompassed -5.0< 77 <-3.0 and 3.5< 77 <5.25, respec- 
tively. In the mid-77 range, the regions covered by the 
OCT subdetector(OCT- and OCT+) were used. For the 
determination of v 2 in the positive (negative) 77 region of 
the detector, OCT- (OCT+) was used as the subevent 
region to evaluate ip2 ■ Fiducial cuts in 77 were used to 
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avoid acceptance edge effects on the flow signal. 

Even after resolution correction, Monte Carlo simu- 
lations showed a residual suppression of the flow signal 
from background particles carrying no flow information, 
the so-called non-flow background. This effect was stud- 
ied in detail using simulated data with a known amount of 
flow. By comparing the output resolution corrected flow 
signal to the input flow signal for many samples of simu- 
lated data with different forms and magnitudes of input 
flow (V2), suppression correction factors, C, were deter- 
mined for each bin of centrality and n in the analysis. 
The suppression factors were found to be independent of 
the assumed flow magnitude and its form. Furthermore, 
as a function of 77, the correction is a constant 12% in the 
OCT and ranges from 15-30% in most of the rings, where 
the backgrounds are higher. As a function of centrality 
the correction is a flat 12%. The final corrected value of 
V2 was determined by 



V2(r?, centrality) 



V2 (n, centrality) 



(7) 



This quantity, averaged over the region -1.0< n <1.0, is 
presented in Fig. [3] as a function of the number of par- 
ticipants, N par t. The n dependence of V2, event averaged 
over centrality, is shown in Fig. ||. For the sample of 
accepted events entering this plot, < N part >=191. This 
is substantially more than the overall < N part >=98 for 
minimum bias events in Hijing, reflecting the trigger inef- 
ficiency and the centrality bias in our vertex reconstruc- 
tion. 
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FIG. 4. Elliptic flow, averaged over centrality, as a func- 
tion of 77. The black error bars represent la statistical errors 
and the gray bands represent the systematic error at ~90% 
confidence level. The points on the negative side are reflected 
about r?=0 and shown on the positive side as open circles. 

Numerous sources of systematic error were investi- 
gated, including effects due to the energy cut, hit merg- 
ing, subevent definition, knowledge of the beam orbit 
relative to the detector, shape of the dN/d?y distribu- 
tion, vertexing algorithm, transverse vertex cuts, mag- 
netic field configuration and suppression correction de- 




termination. The effect of these sources depended both 
on n and centrality. In general, the systematic error aris- 
ing from each source was determined by conservatively 
varying that specific aspect of the analysis (or several as- 
pects in concert) and quantifying the change in the final 
V2 result as a function of n and centrality. The individual 
contributions were added in quadrature to derive the 90% 
confidence level error shown as the gray band in Figs. |3| 
and |]. In addition to the procedure described above, 
the systematic error on the lower side of some points in 
Figs. I and I was increased to reflect the reduced sensi- 
tivity of the analysis to the low level of flow in those bins 
as determined from Monte Carlo studies. 

These results represent the first measurement of V2 as 
a broad function of n at RHIC and agree with the mid- 
rapidity measurements of V2 from the STAR (lj] and 
PHENIX |L6| collaborations. The extended angular cov- 
erage portrayed in Fig. || clearly shows a systematic drop 
in the magnitude of v 2 with \n\. This drop could be 
due, in part, to correlations not associated with the reac- 
tion plane, a drop in (pt) due to dynamics or kinematics 
p7| , substantial directed flow, or a change in the under- 
lying mechanism of particle production as a function of 
n. Since these data appeared in a preliminary form [ p8| , 
theoretical efforts have had only limited success in repro- 
ducing the n dependence |l8|,|29|]. These results call into 
question the common assumption of longitudinal boost 
invariance over a large region of rapidity in RHIC colli- 
sions. Further work is needed to develop a clear, three- 
dimensional picture of these collisions. 
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